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I. INTRODUCTION 

Measurements of the parameter sin2/3 [1, 2] have es- 
tablished CP violation in the B° meson system. These 
measurements provide strong support for the Kobayashi 
and Maskawa model of this phenomenon as arising from a 
single phase in the three-generation CKM quark-mixing 
matrix [3] . We present in this paper results from a time- 
dependent analysis of the B° — > 7r + 7r~7r° [4] Dalitz plot 
which is dominated by intermediate vector resonances 
(p) . The goal of this analysis is the simultaneous extrac- 
tion of the strong transition amplitudes and the weak in- 
teraction phase a = arg [— V td V t * b /V ud V* b ] of the Unitar- 
ity Triangle. In the Standard Model, a non-zero value for 
a is responsible for the occurrence of mixing-induced CP 
violation in this decay. The BABAR and Belle experiments 
have obtained constraints on a from the measurement of 
effective quantities sin2a e ff in B decays to tt + tt~ [5, 6] 
and from p + p~ [7, 8], using an isospin analysis [9]. 

Unlike 7r + 7r~, p ± 7r T is not a CP eigenstate and four 
flavor-charge configurations (B°(B Q ) — > /0 ± 7r =F ) must be 
considered. The corresponding isospin analysis [10] is 
unfruitful with the present data sample since two pen- 
tagonal amplitude relations with 12 unknowns have to 
be solved (compared to 6 unknowns for the tt + it~ and 
p + p~ systems). However, it has been pointed out by 
Snyder and Quinn [11] that one can obtain the necessary 
degrees of freedom to constrain a without ambiguity by 
explicitly including in the analysis the variation of the 
strong phases of the interfering p resonances in the Dalitz 
plot. 



DECAY AMPLITUDES 



We consider the decay of a spin-zero B° meson with 
four-momentum ps into the three daughters ir + , ir~ , ir°, 
with and po their corresponding four-momenta. 

We take as the independent (Mandclstam) variables the 
invariant squared masses of the charged and neutral pions 



s+ = (p++pa) 2 



(p-+p ) 2 . (1) 



The invariant squared mass of the positive and negative 
pion, s = (p+ +P-) 2 , is obtained, from energy and 
momentum conservation, 



so = m 2 B o + 2ml + + m^ - s+ - s_ 



(2) 



The differential B° decay rate distribution as a function 
of the variables defined in Eq. (1) (i.e., the Dalitz plot) 



t Also with Universita della Basilicata, Potenza, Italy 

§Also with IPPP, Physics Department, Durham University, 

Durham DH1 3LE, United Kingdom 



reads 

aT{B° -» n+ir-ir ) 



\A 37r \ 



(2tt) 3 32m' 



ds+ds- 



(3) 



where A 3v is the Lorentz-invariant amplitude of the 
three-body decay [12]. 

We assume in the following that the amplitude A 37T and 
its CP conjugate A 3 ^, corresponding to the transitions 
B° — ► 7T + 7r - 7r° and B° — ► 7r + 7r~7r°, respectively, are 
dominated by the three resonances p + , p~ and p . The p 
resonances are assumed to be the sum of the ground state 
p(770) and the radial excitations p(1450) and p(1700), 
with masses and widths determined by a combined fit 
to t+ — > 77 T 7r + 7r° and e+e~ — > tt + tt~ data [13]. Since 
the hadronic environment is different in B decays, we 
do not rely on this result for the relative p(1450) and 
p(1700) amplitudes but instead simultaneously measure 
them with the CP parameters from the fit. Variations 
of the other parameters and possible contributions to 
the B° — > 7T + 7r~7r° decay other than the p resonances 
are studied as part of the systematic uncertainties (Sec- 
tion V). 

We write the A 37T and A 37V amplitudes [11, 14] 



A 3 n = f+A+ + f_A- 
A 3 , = f+A+ + f_A- 



+ foA° 
+ foA° 



(4) 
(5) 



where the f K (with n = {+,—,0} denoting the charge 
of the p from the decay of the B° meson) are functions 
of the Dalitz variables s+ and s_ that incorporate the 
kinematic and dynamical properties of the B° decay into 
a vector p resonance and a pseudoscalar pion. The A K 
are complex amplitudes that include weak and strong 
transition phases and that are independent of the Dalitz 
variables. 

Following Ref. [13], the p resonances are parameterized 
in f K (where k is the charge) as a sum of the p(770), 
p(1450), and p(1700) resonances: 

f K (s) cx F p(770) ( S ) + a p ,e^ P 'F p(U5Q) ( S )+ (6) 
apue^p" F p{l700) {s) , 

where the F p are modified relativistic Breit-Wigner func- 
tions introduced by Gounaris and Sakurai (GS) [15] and 
the a p (4> p ) are the magnitudes (phases) of the higher 
mass p resonances, relative to the p(770). In this anal- 
ysis, we assume that the a and <j) for / + and /_ are the 
same while for /o, corresponding to the much smaller p° 
component, we fix a' p and a" p to zero. Note that the defi- 
nitions (4) and (5) are based on the assumption that the 
relative phases between the p(770) and its radial excita- 
tions are CP-conserving. 

Due to angular momentum conservation, the spin-one 
p resonance is restricted to a helicity-zero state. For a p K 
resonance with charge k, the GS function is multiplied by 
the kinematic function — 4|pJ|p T | cos^, where p K is the 
momentum of either of the daughters of the p resonance, 
p T is the momentum of the particle not from the p decay, 
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and cos9 K is the cosine of the helicity angle of the p K 
all defined, in the p-resonance rest frame. For the p + 
(p~), + (6-) is defined by the angle between the ir° 
(n~) momentum in the p + (p~) rest frame and the p + 
(p~) flight direction in the B° rest frame. For the p°, 9q 
is defined by the angle between the ir + momentum in the 
p° rest frame and the p° flight direction in the B° rest 
frame. With these definitions, each pair of GS functions 
interferes destructively at equal masses. 

The factor of cos8 K in the kinematic functions leads 
to an increased population in the interference regions be- 
tween the different p bands in the Dalitz plot, and thus 
increases the sensitivity of this analysis [11]. 



The 27 coefficients (9)-(14) are real-valued parameters 
that multiply the f K f* bilinears (where k and a denote 
the charge of the p resonances) [16]. These coefficients 
are the observables that are determined by the fit. Each 
of the coefficients is related in a unique way to physically 
more intuitive quantities, such as tree-level and penguin- 
type amplitudes, the angle a, or the quasi-two-body CP 
and dilution parameters [17] (cf. Section VI). The pa- 
rameterization (8) is general: the information on the mir- 
ror solutions (e.g., on the angle a) that are present in the 
transition amplitudes A K , A K is conserved. 

The decay rate (7) is used as a probability density 
function (PDF) in a maximum-likelihood fit and must 
therefore be normalized: 



B. TIME DEPENDENCE 

With At = t 3n —t t3xg defined as the proper time interval 
between the decay of the fully reconstructed B 3n and 
that of the other meson _B° ag from the T(AS), the time- 
dependent decay rate |^^.(At)| 2 (1^4.^ (A.£) | 2 ) when the 
B° is a B° (B°) is given by 



e -\At\/r B0 

4t b o 
T(|A 37r | 2 - 

±2Im 



^37r| 2 + |^3n-| 2 

|I 37r | 2 ) cos(Am d Ai) 
sm(AnidAt) 



-An A* 



(7) 



where r^o is the mean neutral B lifetime and Arrid is 
the B°B° mass difference. Here, we have assumed that 
CP violation in B°B° mixing is absent (\q/p\ = 1) and 
the lifetime difference between B^ and Bl is Ar^ = 0. 
Inserting the amplitudes (4) and (5), one obtains for the 
terms in Eq. (7) 

|A 3 .| 2 ±L4 3 .| 2 = £ \f K \ 2 Ut + 

«e{+,-,o} 

E 2 ( Re [/„/;] Ut a Kc Im [/„£] U^) , 

K<cr£{ + ,-,0} 



Im 



E i/«i 2/ « -» 

«£{ + ,-, 0} 

E (Re[f K r a }i^+im[u:]i^) 

K<ae{ + ,-,0} 



(8) 



with 





= \A K \ 2 ±\A K 


I 2 , 


(9) 


rr±,Re 

^ K<J 


= Rc[A K A a * 


±A K A°*] , 


(10) 


jj±,Im 


= lm[A K A cr * 


±A K A a *] , 


(11) 




= lm[A K A K * 




(12) 


tRc 

K<J 


= Rc^A 17 * 


-A" A"*] , 


(13) 


rim 

KG 


= Im [I K A" 


+ A a A K *] . 


(14) 



where 



+ \A 3 ^) 



\AUAt)\\ (15) 



|^| 2 + |^3.| 2 ) 



E <i/«i 2 >^ + 

k£{+,-,0} 

+ 2Re E (/«/:) + * ' ^ Im ) , (16) 

k<ct£{ + ,-,0} 

where (...) denotes the expectation value over the Dalitz 
plot. The complex expectation values (f K f£) are ob- 
tained from Monte Carlo integration of the Dalitz 
plot (3), taking into account acceptance and resolution 
effects. In this paper, we determine the relative values of 
U and / coefficients to {/+ leaving 26 free coefficients. 

The choice to fit for the U and / coefficients rather 
than fitting for the complex transition amplitudes and 
the weak phase a directly is motivated by the following 
technical simplifications: (i) in contrast to the ampli- 
tudes, there is a unique solution for the U and I coef- 
ficients requiring only a single fit to the selected data 
sample; (ii) in the presence of background, we find that 
the errors on the U and I coefficients are approximately 
Gaussian, which in general is not the case for the ampli- 
tudes; and (in) the propagation of systematic uncertain- 
ties and the averaging between different measurements 
are straightforward for the U and I coefficients. 

The {/+ coefficients are related to resonance branch- 
ing fractions and charge asymmetries; the U~ coeffi- 
cients determine the relative abundance of the B° decay 
into p + ir~ and p~n + and the time-dependent direct CP 
asymmetries. The I K measure mixing-induced CP viola- 
tion and are sensitive to strong phase shifts. Finally, the 

uta Re ^ Im ^ and I^ Im ^ coefficients describe the interfer- 
ence pattern in the Dalitz plot, and their presence dis- 
tinguishes this analysis from the quasi-two-body analysis 
previously reported in Ref. [17]. They represent the addi- 
tional degrees of freedom that allow one to determine the 
unknown penguin contribution and the relative strong 
phases. However, because the overlap regions of the res- 
onances are small and because the events reconstructed 
in these regions suffer from large misreconstruction rates 
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and background, a substantial data sample is needed to 
perform a fit that constrains all amplitude parameters. 

We determine the physically relevant quantities in a 
subsequent least-squares fit to the measured U and / 
coefficients. 



C. THE SQUARE DALITZ PLOT 

Both the signal events and the combinatorial e + e~ — > 
qq (q = u, d, s, c) continuum background events populate 
the kinematic boundaries of the Dalitz plot due to the low 
final state masses compared with the B° mass. We find 
the representation (3) is inconvenient when one wants to 
use empirical reference shapes in a maximum-likelihood 
fit. Large variations occurring in small areas of the Dalitz 
plot are very difficult to describe in detail. These regions 
are particularly important since it is here that the inter- 
ference between p-resonances, and hence our ability to 
determine the strong phases, occurs. We therefore apply 
the transformation 



ds + ds_ 



I det J I dm' dd' 



(17) 



which defines the Square Dalitz plot (SDP). The new co- 
ordinates are 



1 

— arccos 



rriQ — to™ 111 



- 1 



= -e , (is) 

7T 



where too = ^/sq is the invariant mass of the charged 
particles, TO nax = mgo — m n o and TOq 111 " = 2m^+ are 
the kinematic limits of too, 9q is the p° helicity angle, 
and J is the Jacobian of the transformation. Both new 
variables range between and 1 . The determinant of the 




Jacobian is given by 

i,j ,< i * n *i dm o <9cos6> 
|dctJ| = 4| P+ ||p |too- — ■ aoi , (19) 



dm' 



89' 



where |p^| = \JE* + - m\+ and |pg| = ^ E$ - to^, and 

where the energies E+ and Eq are defined in the ir + ir~ 
rest frame. Figure 1 shows the determinant of the Ja- 
cobian as a function of the SDP parameters ml and 9'. 
If the events in the nominal Dalitz plot were distributed 
according to a uniform (non-resonant) three-body phases 
space, their distribution in the SDP would match the plot 
of | det J|. 

The effect of the transformation (17) is illustrated in 
Fig. 2, which displays the nominal and square Dalitz plots 
for simulated signal events generated with Monte Carlo. 
As is shown, the transformation is benificial because: (i) 
it expands the regions of interference so that equal size 
bins cover this region in more detail; and (n) it avoids 
the curved edge of bins on the boundary. This simulation 
does not take into account any detector effects and cor- 
responds to a particular choice of the decay amplitudes 
for which destructive interferences occur where the p res- 
onances overlap. To simplify the comparison, hatched 
areas showing the interference regions between p bands 
and dashed isocontours ^/s+,-,0 = 1-5 GeV/c 2 have been 
superimposed on both Dalitz plots. 



II. THE BABAR DETECTOR AND DATASET 

The data used in this analysis were collected with the 
BABAR detector at the PEP-II asymmetric-energy e + e~ 
storage ring at SLAC between October 1999 and August 
2006. The sample consists of about 346 fb _1 , correspond- 
ing to (375 ± 4) x 10 6 BB pairs collected at the T(45) 
resonance ( "on- resonance" ) , and an integrated luminos- 
ity of 21.6 fb" 1 collected about 40 MeV below the T(45) 
( "off- resonance" ) . In addition, we use GEANT4 [18] sim- 
ulated Monte Carlo (MC) events to study detector effi- 
ciency and backgrounds. 

A detailed description of the BABAR detector is pre- 
sented in Ref. [19]. The tracking system used for charged 
particle and vertex reconstruction has two main compo- 
nents: a silicon vertex tracker (SVT) and a drift cham- 
ber (DCH), both operating within a 1.5-T magnetic field 
generated by a superconducting solenoid. Photons are 
identified in an electromagnetic calorimeter (EMC) sur- 
rounding a detector of internally reflected Cherenkov 
light (DIRC), which associates Cherenkov photons with 
tracks for particle identification (PID). Muon candidates 
are identified with the use of the instrumented flux return 
(IFR) of the solenoid. 



FIG. 1: Jacobian determinant (19) of the transformation (17) 
defining the square Dalitz plot (SDP). Such a distribution 
would be obtained in the SDP if events were uniformly dis- 
tributed over the nominal Dalitz plot. 



III. ANALYSIS METHOD 

The U and / coefficients and the B° — > 7r + 7r _ 7r° event 
yield are determined by a maximum-likelihood fit of the 
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FIG. 2: Nominal (left) and square (right) Dalitz plots for Monte-Carlo generated B a —> tt + it~tt u decays. Comparing the two 
Dalitz plots shows that the transformation (17) indeed homogenizes the distribution of events, which are no longer near the 
plot boundaries but rather cover a larger fraction of the physical region. The decays have been simulated without any detector 
effect and the amplitudes A + , A~ and A have all been chosen equal to 1 in order to have destructive interferences where the 
p bands overlap. The main overlap regions between the p bands are indicated by the hatched areas. Dashed lines in both plots 
correspond to = 1.5 GeV/c 2 ; the central region of the Dalitz plot (defined by requiring that all 3 two-body invariant 

masses exceed this threshold) contains few signal events. 



signal and background model to the selected candidate 
events. Kinematic and event shape variables exploiting 
the characteristic properties of the events are used in the 
fit to discriminate signal from background. 



A. EVENT SELECTION AND BACKGROUND 
SUPPRESSION 

We reconstruct B° — ► ir + ir~ir° candidates from pairs 
of oppositely-charged tracks and a 7r° — > 77 candi- 
date. In order to ensure that all events are within the 
Dalitz plot boundary, we constrain the three-pion invari- 
ant mass to the B mass after final selections have been 
made. The largest source of background is from contin- 
uum e + e~ — * qq production. We use information from 
the tracking system, EMC, and DIRC to remove tracks 
for which the PID is consistent with the electron, kaon, 
or proton hypotheses. In addition, we require that at 
least one track has a signature in the IFR that is incon- 
sistent with the muon hypothesis. This selection retains 
92% of signal events while rejecting 42% of continuum 
background events. The 7r° candidate mass 771(77) must 
satisfy 0.11 < 777(77) < 0.16GeV/c 2 , where each photon, 
7, is required to have an energy greater than 50MeV in 
the laboratory frame (LAB) and to exhibit a lateral pro- 
file of energy deposition in the EMC consistent with an 
electromagnetic shower. 

A i?-meson candidate is characterized kinemati- 



cally by the beam-energy substituted mass toes = 
V^blLJ 2 - (pT) 2 and energy difference AE = E* B - 
\^fs, where (Eb,Pb) and (Eq,Po) are the four- vectors of 
the B-candidate and the initial electron-positron systems 
respectively. The asterisk denotes the center-of-mass 
(CM) frame and s is the square of the CM energy. We re- 
quire 5.272 < toes < 5.288 GeV/c 2 , which retains 81% of 
the signal and 8% of the continuum background events. 
The AE resolution exhibits a dependence on the 7r° en- 
ergy and therefore varies across the Dalitz plot. To avoid 
bias in the Dalitz plot, we introduce the transformed 
quantity AE' = (2AE - AE+ - AE_)/(AE + - AE-), 
with AE±(m ) — c± — (c± =F c) (too/to™ 3 -*) 2 , where 
too = ^/sq is strongly correlated with the energy of the 

= -0.140GeV, 
require —1 < 



We use the values c = 0.045 GeV, c_ 



= 0.080 GeV, 



= 5.0 GeV, and 



AE < 1. These values have been obtained from Monte 
Carlo simulation. The requirement retains 75% (25%) of 
the signal (continuum) events. 

Backgrounds arise primarily from random combina- 
tions of 7t ± and 7T° candidates in continuum events. Con- 
tinuum events tend to have a more "jet-like" structure 
than B decays which are produced nearly at rest in the 
CM system. To enhance discrimination between signal 
and continuum, we use a neural network (NN) [20] to 
combine four discriminating variables: the angles with 
respect to the beam axis of the B momentum and B 
thrust axis in the T(AS) frame, and the zeroth and sec- 
ond order polynomials Lq,2 of the energy flow about 
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the B thrust axis. The polynomials are defined by 
L n = J2 i p i -\cos9 i \ n , where 6^ is the angle with respect to 
the B thrust axis of any track or neutral cluster i, pi is its 
momentum, and the sum excludes the B candidate. The 
NN is trained with off-peak data and simulated signal 
events. The final sample of signal candidates is selected 
with a requirement on the NN output that retains 77% 
(8%) of the signal (continuum) events. A total of 35444 
on-peak data events pass the selection. 

The time difference At is obtained from the measured 
distance between the z positions (along the beam direc- 
tion) of the B®^ and £?° ag decay vertices, and the boost 
/?7 = 0.56 of the e+e" system: At = Az/f3~/c. The 
f?° ag vertex is determined from the charged particles in 
the event not included in the signal B. To determine 
the flavor of the B® ag we use the B flavor-tagging algo- 
rithm of Ref. [1]. This produces six mutually exclusive 
tagging categories. We improve the efficiency of the sig- 
nal selection by retaining untagged events in a seventh 
category which contribute to the measurement of direct 
CP violation. 

Multiple B candidates passing the full selection occur 
in 16% (p ± tt z f) and 9% (p 7r°) of p(770) MC events. If 
the multiple candidates have different tt° candidates, we 
choose the B candidate with the reconstructed tt° mass 
closest to the nominal ir° mass; in the case that more 
than one candidate have the same ir°, we arbitrarily chose 
a reconstructed B candidates passing the selection (this 
occurs in 4% of events). 

The signal efficiency determined from MC simulation 
is 24% for B° -> p ± ^ and B° -> p°7r° events, and 
11% for non-resonant B° — > -K + ir~ir° events. The signal 
efficiency distribution on the SPD is shown in Figure 3. 




0.2 0.4 0.6 0.8 
m' 



FIG. 3: The signal efficiency distribution on the square Dalitz 
plot. Note that the plot is folded in 6' since the distribution 
is nearly symmetric in this variable. 

The signal events passing the event selection are a com- 
bination of correctly reconstructed ( "truth-matched" , 
TM) events and mis-reconstructed ("self-cross-feed", 



SCF) events. Of the selected signal events, 22% of 
B° -» p ± n*, 13% of B° -> p°ir°, and 6% of non- 
resonant events are mis-reconstructed, according to MC. 
Mis-reconstructed events occur when a track or neutral 
cluster from the tagging B is assigned to the recon- 
structed signal candidate. This occurs most often for 
low-momentum particles and photons; hence the mis- 
reconstructed events are concentrated in the corners of 
the standard Dalitz plot. Since these are also the ar- 
eas where the p resonances overlap strongly, it is im- 
portant to model the mis-reconstructed events correctly. 
The details of the model for the distributions of mis- 
reconstructed events in the Dalitz plot are described in 
Section III CI. 



B. BACKGROUND FROM OTHER B DECAYS 

We use MC simulated events to study the background 
from other B decays. More than one-hundred chan- 
nels were considered in these studies, of which 29 are 
included in the final likelihood model. These exclusive 
S-background modes are grouped into eighteen differ- 
ent classes according to their kinematic and topological 
properties: six for charmless B + decays, eight for charm- 
less B° decays and four for exclusive charmed B° decays. 
Two additional classes account for inclusive B° and B + 
charmed decays. 

Table I summarizes the twenty background classes that 
are used in the fit. For each mode, the expected num- 
ber of selected events is computed by multiplying the 
selection efficiency (estimated using MC simulated de- 
cays) by the branching fraction, scaled to the dataset 
luminosity (346 fb _1 ). The world average branching ra- 
tios have been used for the experimentally known decay 
modes[12, 21]. When only upper limits are given, they 
have been translated into branching ratios including ad- 
ditional conservative hypotheses (e.g., 100% longitudinal 
polarization for B — > pp decay) if needed. 



C. THE MAXIMUM-LIKELIHOOD FIT 

We perform an unbinned extended maximum- 
likelihood fit to extract the total B° — > 7r + 7r~7r event 
yield, and the U and I coefficients defined in Eqs. (9)- 
(14). The fit uses the variables At, to', 6' ', Toes, AE' , and 
NN output to discriminate signal from background. The 
At distribution is sensitive to mixing-induced CP viola- 
tion but also provides additional continuum-background 
rejection. 

The selected on-resonance data sample is assumed 
to consist of signal, continuum-background, and B- 
background components, separated by the flavor and tag- 
ging category of the tag side B decay. The probability 
density function for event i in tagging category c is 
the sum of the probability densities of all components, 
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TABLE I: Summary of the B-background modes taken into account for the likelihood model. They have been grouped in 
20 classes: charmless B + (six), charmless B° (eight), exclusive charmed B° (four) and inclusive 5° and charmed B + decays. 
Modes with at least two events expected after final selection have been included. 



Class 


Mode 




Al/A.JJLL- ILU llUlilUtl (Jl C V CIILjO 


1 


7-1-+- -4- n 

B + - pV[long] 


19.1 ± 3.5 


57 ± 11 


1 


-A- / / \ -I- \ O 

B + -» of (-» (pn) + )n" 


20.0 ± 15.0 


35 ± 26 


i 
i 


D+ . ,,0/ . „±„T\„+ 
D — » <2l( — » p 7T T j7T 


on n _i_ 1 k n 
20.0 ± lo.O 


21 ± 15 


2 


D+ + 


8.7 ± 1.0 


81 ± 9 


3 




4.3 ± 0.6 


7± 1 


3 


_D — ► 7T XV 5 ^ — ► 7T 7T J 


O.o zt U.4 


11-1- 1 

11 zt 1 


4 




lU.o zt 1.4 


/u zt y 


4 


_D — > 7T — ' TV TV J 


O.I T U.Z 


17-1- 9 
1 / zt z 


5 


D+ + 


U.tJ zc u.u 


1 ^ 4- 9 
_i_tj zrz Zi 


5 


B + K+7V 


12.1 ± 0.8 


9± 1 


6 


D+ . ( Tfi**) (1 /IQnW"^ k f ^+-77-77- "l + 

_D — ► V-"- ^l^OUJvr^ — > ^A 7T7V ) 


zy.u zt 0.4 


/19 -1- 
4Z Zt D 


7 


JD — > 7V A ^ — > i\ gTT j 


o.o zt U.4 


9-1- 1 
Z zt 1 


00 00 


R0 _^ „+ - 

R° — > (n, it) 
D — > \CLLTl) 


zo.z zh O. ( 

jy. / zt o. ( 


^4 zt 11 

4o Zt 4 


9 


D — > A 7T 


lo.y zh U. ( 


lO zt 1 


10 


i3 — > 7T l\ [ — > A 7V J 


o.o zt U.4 


99 -1- 9 
ZZ zt Z 


10 


D — > A v '(14oU)7T — > A7T7T 


11. z dz z.z 


Z04 it Oi 


11 
11 


d — > 7 A (,o»2, 14oUJ(^ — > (A tt ) ) 
R° — > 7T <T* ('^ <T+7r" % l 


07 /I _|_ 1 C 

Z ^.4 zt 1.0 

i o _i_ r\ r 
l.o zrz u.o 


1 C _|_ 1 

10 zt 1 

1 4- zl 

iu zrz ^ 
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0.4 ± 0.2 


3± 2 


12 




9.9 ± 1.6 


114 ± 17 


13 


B ° -» # +7r ~ 7r °[nonrc S ] 


4.6 ± 4.6 


42 ± 38 


14 




5.8 ± 0.5 


55 ± 4 


15 


B°^D-(^7r-7r°)7r+ 


7.5 ± 2.3 


661 ± 203 


16 


B° - D°(-^ K+n-)n° 


11.0 ± 3.2 


110 ± 32 


17 


B°^D°(-^ 7v+7v-)n" 


0.4 ± 0.1 


39 ± 10 


18 


B° -> J/^H e + e~,p+p-)7r () 


2.6 ± 0.5 


85 ± 17 


19 
20 


B u — > {neutral generic b — > c decays} 
B + — > {charged generic & — > c decays} 




188 ± 30 
431 ± 40 



namely 



-TM,i 



/sCF^STr-SCF^ 



+ ^2 N B°jfl3°j'PB°, 



(20) 



where is the total number of tt + tt it signal events in 
the data sample; is the fraction of signal events that 
are in tagging category c; J SCF is the fraction of SCF 
events in tagging category c, averaged over the Dalitz 
plot; "P^-tm i and V^-scf i are the products of PDFs 
of the discriminating variables used in tagging category c 
for TM and SCF events, respectively; is the number 
of continuum events that are in tagging category c; <7tag,i 



is the tag flavor of the event, defined to be +1 for a 
.B° ag and —1 for a B® ag ; Aq^tag parameterizes possible 
flavor tag asymmetry in continuum events; P%g i is the 

continuum PDF for tagging category c; N^ ss (N^ ss ) is 
the number of charged (neutral) _B-related background 
classes considered in the fit; N B +j (Nsoj) is the number 
of expected events in the charged (neutral) _B-background 
class j; fg+j (/soj) is the fraction of charged (neutral) B- 
background events of class j that are in tagging category 
c; A B + .tag..? describes a possible flavor tag asymmetry in 
the B + background class j; Pg+ 4 • is the £? + -background 
PDF for tagging category c and class j; and Pg ^ is 
the neutral-B-background PDF for tagging category c 
and class j. Correlations between the flavor tag and the 
position in the Dalitz plot are absorbed in tag-flavor- 
dependent Dalitz plot PDFs that are used for B + and 
continuum background. 



The PDFs P c x {X = {TM,SCF,qq,B+/B }) are the 
product of the four PDFs of the discriminating variables, 
x\ = niESi x 2 = Ai?', and x% = NN output, and the 
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triplet x A = {m',d',At}: 

fe=i 

where i is the event index and j is a _B-background class. 
The extended likelihood over all tagging categories is 
given by 

7 N c 

£^U eN Il V *> ( 22 ) 

c=l i=l 

where A^ c is the total number of events expected in cat- 
egory c. 

A total of 68 parameters, including the inclusive signal 
yield N 37T and the 26 U and I coefficients from Eq. (7), 
are varied in the fit. Most of the parameters describing 
the continuum distributions are also free in the fit. The 
parameterizations of the PDFs are described below and 
are summarized in Tab. II. 



1. THE At AND DALITZ PLOT PDFS 

The Dalitz plot PDFs require as input the Dalitz plot- 
dependent relative selection efficiency e — s(m',6'), and 
the SCF fraction, /scf = /scf(to', 9'). Both quantities 
are taken from MC simulation. Away from the Dalitz 
plot corners the efficiency is uniform, while it decreases 
when approaching the corners where one of the three 
particles in the final state is almost at rest in the LAB 
frame so that the acceptance requirements on the particle 
reconstruction become restrictive. Combinatorial back- 
grounds, and hence SCF fractions, are large in the cor- 
ners of the Dalitz plot due to the presence of soft neutral 
clusters and tracks. 

For an event i, we define the time-dependent Dalitz 
plot PDFs 

i$r-TM,i = £i (1 - /scf,*) I det J *\ l^(Ai)| 2 ,(23) 

^r-scF.i = d /scf,* I det J *l l^(Ai)| 2 , (24) 

where Pz-n-TM.i and Pstt-scf,* are normalized. The 
normalization involves the expectation values (s (1 — 
/sCF)|detJ|.r.r*) and (e/sc F |detJ|/ K .r*) for TM 
and SCF events, where the indices n, a run over all res- 
onances belonging to the signal model. The expectation 
values arc modcl-dcpcndcnt and are computed with the 
use of MC integration over the square Dalitz plot: 

( £ (i-/sc F )idetj|rr*) = 

Jo Jq s (1 - / SCF ) I det J\ rr* dm'dff 
Jo Jo s I det J| dm'dff 

and similarly for (/scf£ | det J| f K f (T *), where all quan- 
tities in the integrands are Dalitz-plot dependent. 



Equation (20) invokes the phase space-averaged SCF 
fraction / SCF = (/ SCF | det J\ f K f a *). The PDF normal- 
ization is decay-dynamics-dependent and is computed it- 
eratively. We determine the average SCF fractions sepa- 
rately for each tagging category from MC simulation. 

The width of the dominant p(770) resonance is large 
compared to the mass resolution for TM events (about 
8MeV/c 2 Gaussian resolution). We therefore neglect 
resolution effects in the TM model. Mis-reconstructed 
events have a poor mass resolution that strongly varies 
across the Dalitz plot. These events are described in the 
fit by a two-dimensional resolution function 

R SCF (m' r ,e' r ,m' t ,e' t ) , (26) 

which represents the probability to reconstruct at the 
coordinate (m' r , 9' r ) an event that has the true coordinate 
(m' t ,9' t ). This function obeys the unitary condition 

l l 

J J RscF(m' r , 9' r , m' t , 9' t ) dm' r d9' r = 1, (27) 
o o 

and is convolved with the signal model. The -Rscf func- 
tion is obtained from MC simulation. 

The dynamical information in the signal model is de- 
scribed in Section I A and is connected with At via the 
matrix element in Eq. (7), which serves as the PDF. The 
PDF is modified by the effects of mistagging and the lim- 
ited vertex resolution [17]. The At resolution function for 
signal and i?-background events is a sum of three Gaus- 
sian distributions, with parameters determined by a fit 
to fully reconstructed B° decays [1]. Since the major- 
ity of SCF events arise from mis-reconsructed 7r° decays 
which do not affect the vertex resolution, we use the same 
resolution function for TM and SCF events. 

The Dalitz plot- and Ai-dependent PDFs factorize for 
the charged- B background modes, but not necessarily for 
the B° background due to B°B° mixing. 

The charged £?-background contribution to the likeli- 
hood (20) involves the parameter A B +^ tag , multiplied by 
the tag flavor q tag of the event. In the presence of sig- 
nificant "tag-'charge" correlation (represented by an ef- 
fective flavor tag versus Dalitz coordinate correlation), it 
parameterizes possible fake direct CP violation or asym- 
metries due to detector effects in these events. We o also 
use separate square Dalitz plot PDFs for B° and B fla- 
vor tags, and a flavor-tag-averaged PDF for untagged 
events. The PDFs are obtained from MC simulation and 
are described with the use of non-parametric functions. 
The At resolution parameters are determined by a fit to 
fully reconstructed B + decays. For each £? + -background 
class we obtain effective lifetimes from MC to account 
for the mis-reconstruction of the event that modifies the 
nominal At resolution function. 

The neutral-i? background is parameterized with 
PDFs that depend on the flavor tag of the event. In the 
case of CP eigenstates, correlations between the flavor 
tag and the Dalitz coordinate are expected to be small. 
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TABLE II: Summary of PDF parameterizations where G=Gaussian, PX 
biCB=bifurcated Crystal Ball. See Section III C 1 for a detailed description 


=X-order polynomial, NP=non-parametric, and 
af the Dalitz plot parameterization for signal. 


Variable 


TM Signal 


SCF Signal 


Continuum 


B-Background 


AE 

tubs 

Neural Net 

Dalitz 

At 


GG 
biCB 

NP 
see text 
GGG 


G 
NP 
NP 

see text 
GGG 


P2 

Argus 
P3 
NP 

GGG 


NP 
NP 

NP 
NP 

GGG 



However, non-CP eigenstates, such as afn T , may exhibit 
such correlations. Both types of decays can have direct 
and mixing-induced CP violation. A third type of de- 
cay involves charged kaons (e.g. p ± K T ) and docs not 
exhibit mixing-induced CP violation, but usually has a 
strong correlation between the flavor tag and the Dalitz 
plot coordinate, because these decays correspond to B- 
flavor eigenstates. The Dalitz plot PDFs are obtained 
from MC simulation and are described with the use of 
non-parametric functions. For neutral-i? background, 
the signal At resolution model is assumed. 

The Dalitz plot treatment of the continuum events is 
similar to that used for charged-i? background. The 
square Dalitz plot PDF for continuum background is 
obtained from on-resonance events selected in the Toes 
sidebands (defined as 5.225 < toes < 5.265) and cor- 
rected for a 5% feed-through from B decays. A large 
number of cross checks have been performed to ensure 
the high fidelity of the empirical shape parameterization. 
The continuum At distribution is parameterized as the 
sum of three Gaussian distributions with common mean 
and three distinct widths. The widths scale with the es- 
timated At uncertainty for each event. This yields six 
shape parameters that are determined by the fit. The 
model is motivated by the observation that the At av- 
erage is independent of its error, and that the At RMS 
depends linearly on the At error. 



2. PARAMETERIZATION OF THE OTHER 
VARIABLES 

The toes distribution of TM signal events is parame- 
terized by a bifurcated Crystal Ball function [22], which 
is a combination of a one-sided Gaussian and a Crystal 
Ball function, given as: 



m = { 



' Ce (x-m) 2 /2s 2 R 
Ce (x-rn) 2 /2s 2 L 



for (x — to) > 0, 
for > 2^ > -A, 



pend linearly on to-q. The parameters of the narrow Gaus- 
sian are determined by the fit to data while the others 
are obtained from signal MC. Mis-reconstructed events 
are parameterized by a broad single Gaussian function 
whose parameters are taken from signal MC. 

Both toes an d AE PDFs are parameterized by non- 
parametric functions for all £?-background classes. Con- 
tinuum events are parameterized with an Argus shape 
function [24] 



(28) 



C&e'^ (^- A - S 7T L ) for 3^* < -A. 



The peak position of this function, to, is determined by 
the fit to on-peak data while the other parameters are 
taken from signal MC. A non-parametric function [23] is 
used to describe the SCF signal component. 

The AE' distribution of TM events is parameterized by 
a double Gaussian function, where all five parameters de- 



/(m ES ) = C 



TOES 



'1- 



III 



ES 



TOES 



TO 



ES 



-£(1- 



TO-ES 

™ES 



(29) 

and a second-order polynomial in AE' , with parame- 
ters determined by the fit. The value of m^f is 5.2886 
GeV/c 2 . 

We use non-parametric functions to empirically de- 
scribe the distributions of the NN outputs found in the 
MC simulation for TM and SCF signal events, and for 
B-background events. We distinguish tagging categories 
for TM signal events to account for differences observed 
in the shapes. 

The continuum NN distribution is parameterized by 
a third-order polynomial. The coefficients of the polyno- 
mial are determined by the fit. Continuum events exhibit 
a correlation between the Dalitz plot coordinate and the 
inputs to the NN. To account for this correlation, we 
introduce a linear dependence of the polynomial coeffi- 
cients on the distance of the Dalitz plot coordinate from 
kinematic boundaries of the Dalitz plot. The parameters 
describing this dependence are determined by the fit. 



IV. FIT RESULTS 



event yield of = 2067 ± 86, where the error is sta- 
tistical only. The results for the U and I coefficients are 
given together with their statistical and systematic er- 
rors in Table III. The corresponding correlation matrix 
is given in Table IV. We have generated a sample of 
Monte Carlo experiments to determine the probability 
density distributions of the fit parameters. Within the 
statistical uncertainties of this sample we find Gaussian 
distributions for the fitted U and I coefficients. This al- 
lows us to use the least-squares method to derive other 
quantities from these (Section VI). 
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FIG. 4: Distributions of (a-f) AE' , 771es, NN output, At /a (At), m' , and 0' for samples enhanced in B° — > 7r + 7r _ 7r° signal. The 
dots with error-bars correspond to the on-resonance data. The solid histogram shows the projection of the fit result. The dark 
and light shaded areas represent the contribution from continuum and B-background events respectively. The misreconstructed 
signal events are represented by the dashed histogram. The ratios of signal events over background events in the sample differs 
for each plot, but are typically ~ 1.2. 
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TABLE III: Fit results for the U and 7 coefficients. The errors given are statistical (first) and systematic (second). The free 

normalization parameter is fixed to 1. The coefficients are defined in Eqn. 8. 

Parameter Description Result 

U+ Coefficient of | /+ 1 2 1.0 (fixed) 

C7 + Coefficient of j/„| 2 0.28 ± 0.07 ± 0.04 

UZ Coefficient of | /_ | 2 1.32 ± 0.12 ± 0.05 

Uq Coefficient of |/ | 2 cos(Am d At) -0.03 ± 0.11 ± 0.09 

UZ Coefficient of |/_ | 2 cos(Am d Ai) -0.32 ± 0.14 ± 0.05 

Ul Coefficient of |/+| 2 cos(Am d At) 0.54 ± 0.15 ± 0.05 

7 Coefficient of |/ | 2 sin(Am d At) 0.01 ± 0.06 ± 0.01 

I- Coefficient of j/_j 2 sin(Arn d Ai) -0.01 ± 0.10 ± 0.02 

7+ Coefficient of |/+| 2 sin(Am d Ai) -0.02 ± 0.10 ± 0.03 

U+'l™ Coefficient of Im [/+/*] -0.07 ± 0.71 ± 0.73 

U^ Rc Coefficient of Re [/+/*] 0.17 ± 0.49 ± 0.31 

UZ'l™ Coefficient of Im [/+/*] cos(Am d Ai) -0.38 ± 1.06 ± 0.36 

UZ 3c Coefficient of Re[/+/*] cos(Am d At) 2.23 ± 1.00 ± 0.43 

il^L Coefficient of Im [/+/*] sm(Am d At) -1.99 ± 1.25 ± 0.34 

j£l Coefficient of Re[/+/*] sin(Am d Ai) 1.90 ± 2.03 ± 0.65 

t/+J Im Coefficient of Im[/+/o] -0.16 ± 0.57 ± 0.14 

C/+f c Coefficient of Re[/+/ * ] - 1.08 ± 0.48 ± 0.20 

(7+o Im Coefficient of Im[/+/ *] cos(Am d At) -1.66 ± 0.94 ± 0.25 

UIa Rc Coefficient of Re[/+/ * ] cos(Am d Ai) -0.18 ± 0.88 ± 0.35 

Coefficient of Im [/+/£] sin(Am d At) -0.21 ± 1.06 ± 0.25 

Coefficient of Re[/+/ *] sin(Am d Ai) 0.41 ± 1.30 ± 0.41 



'+o 

,Im 



UZo Coefficient of Im[/_/o] -0.17 ± 0.50 ± 0.23 

UZo 1 " Coefficient of Re[/_ f,*] -0.36 ± 0.38 ± 0.08 

[C(, Im Coefficient of Im[/_/ *] cos(Am d At) 0.12 ± 0.75 ± 0.22 

UZi He Coefficient of Re[/_/ *] cos(Am d Ai) -0.63 ± 0.72 ± 0.32 

755 Coefficient of Im[/_/ *] sin(Am d At) 1.23 ± 1.07 ± 0.29 

7 R ( °, Coefficient of Re[/_/ *] sin(Am d Ai) 0.41 ± 1.30 ± 0.21 
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TABLE IV: Correlation matrix of statistical uncertainties for the U and / coefficients. Since the matrix is symmetric, all 
elements above the diagonal are omitted. 
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The signal is dominated by B° — > p ± n T decays. We 
observe an excess of p°ir° events (see, mainly, Uq), which 
is in agreement with our previous upper limit [25] and 
the latest measurement from the Belle collaboration [26] . 
Wc find the ratio of p(1450)/p(770) ( p(1700)//)(770) ) 
rates to be 0.13 ± 0.04 (0.07 ± 0.04) where the errors 
are statistically only. For the relative strong phase be- 



tween the p(770) and p(1450) (p(1700)) amplitudes we 
find (163±22)° ((5±36)° (statistical errors only), which is 
compatible with the result from r and e + e~ data. These 
results for the p amplitudes are compatible with the find- 
ings in r and e + e~ decays [13]. 

Figure 4 shows distributions of AE', rngs, the NN out- 
put, and At/<r(At), where a{At) is the per-event error 
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FIG. 5: Distribution of the minimum of the three di-pion 
invariant masses, for samples enhanced in B° — > 7r + 7r - 7r° sig- 
nal. The dots with error-bars correspond to the on-resonance 
data. The solid histogram shows the projection of the fit 
result. The dark and light shaded areas represent the con- 
tribution from continuum and B-background events respec- 
tively. The misreconstructed signal events are represented by 
the dashed histogram. The peaks near 0.5 and 1.8GeV/c 2 are 
from B° — ► K^tt and B° — > D + ty~ decays, respectively. 

on At, as well as the Dalitz plot variables m! and 6'. All 
distributions are enhanced in signal content by selecting 
on the ratio of the probability the event is signal to the 
total, P s ig/J2 Pi excluding the variable plotted. Figure 5 
shows the distribution of the minimum of the three di- 
pion invariant masses, again enhanced in signal content. 
This plot shows clearly that p(770) dominates the signal 
component. 

V. SYSTEMATIC STUDIES 

The contributions to the systematic error on the signal 
parameters are summarized in Table V. Table VI sum- 
marizes the correlation coefficients extracted from the 
systematic covariance matrix. For a given systematic 
effect, we vary a parameter in the fit (e.g. the p(77Q) 
mass), refit the data, and construct the systematic co- 
variance matrix for that source based on the deviations 
of the U and / coefficients from the nominal values. The 
matrix element is given as 

Si,j = SiSj (30) 

(31) 

where Si is the difference between the two two fits for 
variable i. The total systematic covariance matrix is 
obtained by adding together the covariance matrices in 
quadrature from the different systematic sources. 

To estimate the contribution to B° — > ir + Tr~ir° decay 



from other resonances and non-resonant decays, we fit 
the on-peak data including these other possible decays in 
the fit model. For simplicity, we assume a uniform Dalitz 
distribution for the non-resonant events and consider pos- 
sible non-p resonances including / (980), /2(1270), and 
a low mass 5-wave er whose mass and width we take to 
be 478MeV/c 2 and 324MeV/c 2 , respectively [12]. The fit 
does not find a significant signal for any of those decays. 
However, the inclusion of the broad, low mass it + it~~ 
S'-wave component significantly degrades our ability to 
identify p°7r° events. The systematic effect (contained in 
the "Dalitz plot model" rows in Table V) is estimated 
by generating Monte Carlo samples including the other 
B° — > 7r + 7T~7r° modes and fitting with the nominal setup, 
where only p(770) is taken into account. 

We vary the mass and width of the p(770), p(1450), 
and /o(1700) resonances within ranges that exceed twice 
the errors found for these parameters extracted from r 
decays and e + e~ annihilations [13], and assign the ob- 
served shifts in the measured U and / coefficients as sys- 
tematic uncertainties ("p,p',p" lineshapc" in Table V). 
Since some of the U and I coefficients exhibit significant 
dependence on the p(1450) and p(1700) contributions, we 
leave their amplitudes (phases and fractions) free to vary 
in all fits. 

To validate the fitting tool, we perform fits on large MC 
samples with the measured proportions of signal, contin- 
uum, and i?-background events. No significant biases are 
observed in these fits. The statistical uncertainties on 
their fit parameters are taken as systematic uncertainties 
("Fit bias" in Table V). 

Another potentially large source of systematic uncer- 
tainty is the i3-background model. The expected event 
yields from the background modes are varied accord- 
ing to the uncertainties in the measured or estimated 
branching fractions ( "-^Background" in Table V). Since 
JS-background modes may exhibit CP violation, the cor- 
responding parameters are varied either within their mea- 
sured ranges (if available) or within ±0.5 (if unmeasured) 
("B background CP" in Table V). 

Other systematic effects are much less important to the 
measurements of U and / coefficients and are combined 
in the "Others" field in Table V. Details are given below. 

The parameters for the continuum events are deter- 
mined by the fit. No additional systematic uncertainties 
are assigned to them. An exception to this is the Dalitz 
plot PDF; to estimate the systematic uncertainty from 
the toes sideband extrapolation, we select large samples 
of off-resonance data by loosening the requirements on 
AE and the NN. We compare the distributions of to' 
and 9' between the Toes sideband and the signal region. 
No significant differences are found. We assign as system- 
atic error the effect seen when weighting the continuum 
Dalitz plot PDF by the ratio of the 2-dimensioal his- 
tograms taken from the signal region and sideband data 
sets. This effect is mostly statistical in origin. 

The uncertainties associated with Arrid and r are esti- 
mated by varying these parameters within the uncertain- 
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TABLE V: Summary of systematic uncertainties. 
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0.047 


0.085 


0.038 


0.062 


0.527 


0.094 


0.039 



ties on the world averages [12]. 

The systematic effects due to the signal PDFs com- 
prise uncertainties in the PDF parameterization, the 
treatment of misreconstructed decays, the tagging per- 
formance, and the modeling of the signal contributions. 

When the signal PDFs are determined from fits to a 
control sample of fully reconstructed B decays to exclu- 
sive final states with charm, the uncertainties are ob- 
tained by varying the parameters within the statistical 



uncertainties. In other cases, the dominant parameters 
have been left free to vary in the fit, and the differences 
observed in these fits are taken as systematic errors. 

The average fraction of misreconstructed signal events 
predicted by the MC simulation has been verified with 
fully reconstructed B — > Dp events [17]. No significant 
differences between data and the simulation are found. 
We vary /scf f° r au tagging categories relatively by 25% 
to estimate the systematic uncertainty. 



19 



TABLE VI: Correlation matrix of systematic uncertainties for the U and I coefficients. Since the matrix is symmetric, all 
elements above the diagonal are omitted. 
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As is done for the signal PDFs, we vary the At resolu- 
tion parameters and the flavor-tagging parameters within 
their uncertainties and assign the differences observed in 
the data fit with respect to the nominal fit as systematic 
errors. 

The systematic errors for the parameters that measure 
interference effects are dominated by the uncertainty in 



the signal model, mainly the description of the p reso- 
nance tails. For the other parameters, the uncertainty 
on the fit bias and the B-background contamination arc 
important. 

As a validation of our treatment of the time depen- 
dence we allow r B a to vary in the fit. We find t b o = 
(1.513 ± 0.066) ps, while the remaining free parameters 
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are consistent with the nominal fit. To validate the SCF 
modeling, we leave the average SCF fractions per tagging 
category free to vary in the fit and find results that are 
consistent with the MC prediction. 



VI. INTERPRETATION OF THE RESULTS 

We can use the results of this time-dependent Dalitz 
analysis to extract the B°(B°) — ► p ± tt t parameters de- 
fined in Ref. [17]: 



± T e -|At|/r 

r Q Z (At) = (i±^-)- ir - 



(32) 



1 + <2ta g (S ± AS) sin(Am d At) 
- Qta g (C ± AC) cos(Am d Ai) 



where Qtag = 1(— 1) when the tagging meson B^ ag is a 
B°(B°). The time- and flavor-integrated charge asym- 
metry A p7T measures direct CP violation and the quanti- 
ties S and C parameterize mixing-induced CP violation 
related to the angle a, and flavor-dependent direct CP 
violation, respectively. The parameters AC p7r and AS pv 
are insensitive to CP violation. 

The U and I coefficients are related to the parameters 
as follows: 



A P ir 







C+ 








S+ 


21+ 
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-uz 




+ UZ 



Cr 



s- 



uz 

2 7_ 



(33) 



where C = (C+ + C~)/2, AC = (C+ - C-)/2, S = 
(S + + S-)/2, and AS = (S+-S-)/2 . The definitions of 
Eq. (33) explicitly account for the presence of interference 
effects, and are thus exact even for a p with finite width, 
as long as the U and I coefficients are obtained with 
a Dalitz plot analysis. This treatment leads to slightly 
increased statistical uncertainties compared to the results 
obtained neglecting the interference effects. 

Using a least-squares method including statistical and 
systematic correlations for the U and / coefficients, we 
obtain: 

A p7r = -0.14 ± 0.05 ± 0.02 , (34) 
C = 0.15 ±0.09 ±0.05 , (35) 
S = -0.03 ±0.11 ±0.04 , (36) 

where the first errors are statistical and the second are 
the systematic uncertainties. For the other parameters 
in the description of the B a (B°) — ► pit decay-time de- 
pendence, we measure 

AC = 0.39 ±0.09 ± 0.09 , (37) 
AS = -0.01 ±0.14 ± 0.06 . (38) 
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FIG. 6: Confidence level contours for the direct CP asymme- 
tries Ap~^ versus ApJ~. The shaded areas represent la, 2a 
and 3a contours, respectively. 



In addition, we measure the B° — > 
parameters and decay fraction to be 



P°tt CP- 



iolati' 



Co 



Si) 



EL 

U + 



/on = 



U+ + UZ + u+ 



-0.10 ±0.40 ± 0.53 , (39) 
0.04 ±0.44 ±0.18 , (40) 
0.136 ±0.036 ±0.039 (41) 



The systematic errors are dominated by the uncertainty 
on the CP content of the £?-related backgrounds. Other 
contributions are the signal description in the likeli- 
hood model (including the limit on non-resonant B — > 
7r + 7T~7r° events), and the fit bias uncertainty. The large 
systematic error on Coo is due to the possible tt + tt~ S- 
wave contribution. The correlation matrix, including sta- 
tistical and systematic uncertainties, of the eight quasi- 
two-body parameters is given in Table VII. 

One can transform the experimentally convenient (un- 
corrected) direct-CP violation parameters C and A p7r 
into Ap~~ , A~+ , defined by 



A + ~ - 



K 



— 12 



a: 



where 



K+-| a + l 
-+12_ 1 



-+|2 



- 1 _ A p7r + C + A p7T AC 
1 + AC + A p7r C 

1 = 1 - AC - Ap*C ' 



(42) 



" = (q/p)A~/A + and k~+ = (q/p)A+/A~, so 
(A~k ) involves only diagrams where the p (it) 
meson is formed from the W boson. We find 



that A+~ 
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TABLE VII: Correlation matrix of the quasi-two-body parameters. 
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A+~ = 0.03 ± 0.07 ± 0.04 , (43) 

A~+ = -0.37 ± 0.16 tHo , (44) 

with a correlation coefficient of 0.62 between A^~ and 
A~+. The confidence level contours including system- 
atic uncertainties are shown in Fig. 6. The significance, 
including systematic uncertainties and calculated by us- 
ing a minimum x 2 method, of direct CP violation is less 
than 3.0ct. 

The measurement of the resonance interference terms 
allows us to determine the relative phase 

(5+_ = arg(A+*A-) (45) 

between the amplitudes of the decays B° — ► p~n + and 
B° — > p + 7r~. Through the definitions in Eqs. (9)- 
(14), we can derive a constraint on <5_| from the mea- 
sured U and I coefficients by performing a least-squares 
minimization with the six complex amplitudes as free 
parameters. The constraint can be improved with the 
use of strong isospin symmetry. The amplitudes A K 
represent the sum of tree-level and penguin-type ampli- 
tudes, which have different CKM factors: the tree-level 
(T K ) B° —> p K ir K transition amplitude is proportional to 
V u dV* b , while the corresponding penguin- type amplitude 
(P K ) involves V q dV* b , where q = u,c,t. Here we denote 
by k the charge conjugate of k, where 7t = when k = 0. 
Using the unitarity of the CKM matrix one can reorga- 
nize the amplitudes and obtain [14] 

A K = T K e - ia + P K , 
(q/p)A K = T*e +ia +P«, (46) 

where the magnitudes of the CKM factors have been ab- 
sorbed in T K , P K , T~ and P K . The Eqs. (46) repre- 
sent 13 unknowns of which two can be fixed due to an 
arbitrary global phase and the normalization condition 
= 1. Using strong isospin symmetry one can iden- 
tify P° = — (P + + P~)/2, which reduces the number of 
unknowns to be determined by the fit to nine. This set 

of parameters provides the constraint on 5-\ , shown in 

the left plot of Fig. 7. We find for the solution that is 
favored by the fit 

<5+_ = (37 ± 37)° , (47) 



where the errors include both statistical and systematic 
effects. There is only a marginal constraint on 5^ ob- 
tained at 95% confidence level (C.L.). 

Finally, following the same procedure, we can also de- 
rive a constraint on a from the measured U and / co- 
efficients. The resulting C.L. function versus a is given 
in the right-hand plot of Fig. 7, including systematic un- 
certainties. Ignoring the mirror solution at a + 180°, we 
find at 68% C.L. 

a = (871S)° • (48) 

Almost no constraint on a is achieved at two sigma and 
beyond. 

VII. SUMMARY 

We have presented a measurement of CP-violating 
asymmetries in B° — ► ir + Tr~ir° decays dominated by the 
p resonance. The results are obtained from a data sample 
of 375 million T(4S) — > BB decays. We perform a time- 
dependent Dalitz plot analysis. From the measurement 
of the coefficients of 26 form- factor bilinears we determine 
the three CP- violating and two CP-conserving quasi-two- 
body parameters, and find no evidence of direct CP vi- 
olation. Taking advantage of the interference between 
the p resonances in the Dalitz plot, we derive constraints 
on the relative strong phase between B° decays to p + n~ 
and /9~7T + , and on the angle a of the Unitarity Triangle. 
These measurements are consistent with the results ob- 
tained by Belle [27] as well as with the expectation of a 
SM fit to all constraints on the CKM matrix [28, 29]. 
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